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Abstract. We study the relaxation dynamics of a quantum Ising chain initially prepared in a product of 
canonical states corresponding each to an equilibrium state of part of the chain at a given temperature. 
We focus our attention on the transverse magnetization for which a general expression is given. Explicite 
results are given for the completely factorized initial state, corresponding to a situation where all the spins 
are thermalized independently, and for the two-temperatures initial state, where part of the chain called 
the system is thermalized at a temperature T s and the remaining part is at a temperature Tb. 

PACS. 75.40.Gb Dynamic properties 

05.70.Ln Non-equilibrium and irreversible thermodynamics 



1 Introduction 

Out of equilibrium properties of quantum chains have 
been the purpose of several recent studies (see Ref. pQ 
for a recent review). In particular, after precursor stud- 
ies |2I.'{I lj at the end of the sixties, some renewed in- 
terest has been put on the relaxation properties of free 
fermionic quantum chains initially prepared in nonequilib- 
rium states |5lfil7l8IHI12ll3| . The influence of randomness 
was also considered in this context (T3|. In this contri- 
bution, we present the relaxation of the transverse mag- 
netization of an Ising quantum chain. The system is ini- 
tially prepared in a factorized state p(0) — Y[j Pj where 

the density matrices pj oc e~P jHj are canonical states of 
parts of the chain thermalized each at a given tempera- 
ture (3J . One can look at that state as obtained from 
a chain initially split into several non-interacting terms: 
Ho = ^2jTij, each in contact with a specific heat bath. 
At time t = 0, the interactions between the different parts 
are switched on, the heat baths removed, so that the time 
evolution of the system is generated quantum mechani- 
cally by the total Hamiltonian H — Ho + J2j Hj,j+i> sucn 
that the system will relax from the initial factorized state 
p(0) toward a new state. 

In this paper, we present the transverse magnetiza- 
tion profiles obtained on two different situations. The first 
considered is a situation where the chain is split into two 
parts, one at low temperature called the "system", the 
other at high temperature called the "bath". So, in the 
system-bath case, we concentrate on the front dynamics 
generated by the interface between the system and the 
bath parts |12I6I13| . The other situation we have consid- 
ered in this study is the droplet-like case, where a low 



temperature finite part of size N is put in contact at both 
ends with high temperature infinite chains. In this case, 
we concentrate on the large time t 3> N properties of 
the transverse magnetization, which shows a scaling be- 
haviour. The paper is organized as follows: in the next sec- 
tion we present the dynamics of the Ising quantum chain 
solved in terms of Clifford operators Then we show 
how we calculate the expectation values of several physi- 
cal quantities, paying a special attention to the transverse 
magnetization. Section 5 contains our results and we end 
with a short summary and discussion of them. 



2 Dynamics 

The Hamiltonian of the system under consideration is 
given by 

k=l k=l 

where the as are the Pauli matrices and where h is the 
transverse field pointing in the z-direction. As it is well 
known, this Hamiltonian can be fermionized thanks to the 
Jordan- Wigner transformation |15) : 



n-l 



3 = 1 

n-l 



so that 



-ir 2 r x 



(2) 

(3) 
(4) 
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where the -Ts are Clifford operators satisfying the anti- 
commutation rules 

{!■;,. I;„\ . (5) 

Introducing the 2L component operator r such that 

rt = (rl,r?,ri,rl,...,rlrl) (6) 

the Hamiltonian is put in the quadratic form 

H = -r ] Tr (7) 

4 

where T is a 2L x 2L matrix. The diagonalization of 
the Hamiltonian can then be performed by the intro- 
duction of the diagonal Clifford operators 7* and j q , re- 
lated to the lattice operators via J 1 * — ^2 q <t>q{n)^ q and 
-Tjj = ipq(n)jq with cj> and ip defined through the eigen- 
value equation 

TV q = e q V q (8) 

with V q (2n — 1) = 4> q (n) and V q (2n) — —iip q (n). It leads 
to the form 

« = «E^. 0) 

The time evolution of the diagonal operators is obtained 
simply by j q (t) = U\{t)^ q U q {t) with 

U q {t) = exp (jf TjT?) = cos f + 7^ sin f . (10) 

Utilising the fact that {7^,7^} = 2Sij5 qq >, we obtain 

2 

^w = Eww>^ ( n ) 

3=1 

where we have defined the pseudo-scalar product as 

(C\D)= l -{C\D) (12) 

with {., .} the anticommutator. The time-dependent lat- 
tice Clifford generators, -T^(i), can then be re-expressed 
in terms of the initial time operators r with the help 
of the inverse transforms 7* = Sfc'WO-^fe aim 7<? = 
J2 k ipq(k)r£. Finally, one obtains 

k.i 

with components 

( r k\ r n(t)) = ^ q (k)<P q (n) cos e q t 

(r^(t)) = (T%\r£(-t)) = -J2Mk)%(n)sme q t 

(r 2 k \rl(t)) = J2Mk)M^)^ S e q t . (u) 



One may notice that the time evolution of the 2L- 
components operator r is given by a 2L-rotation matrix 
R with elements 

JM*) = <W(*)> (15) 

where F^i-x = T} and F^i = rf. The time evolution of r 
is then 

r(t) = R(t)r . (16) 

The set {r£} forms an orthonormal basis of a 2L- 
dimensional linear vector space £ with inner product de- 
fined by (.|.) = h{. J, .}. Hence, every vector X £ E has a 
unique expansion X = k (r k \X)rj:. The string expres- 
sion XiX2-.-X n , with X; € £, is a direct product vector 
of the space £\ (g) £2 (8> ••• <8> £ n which decomposition is 

(17) 



3 Expectation values 

In this work, we have considered systems prepared in fac- 
torized initial states of the form 

p(0) = Y[ Pj (18) 
3 

where the initial density matrix p(0) is given by a product 
of smaller density matrices Pj, corresponding to initially 
non-interacting subsystems. To be more specific, we have 
used the canonical initial distributions 

ft = ^exp(-&W,-) (19) 

where j3j is the inverse temperature of the subsystem Hj 
which is a part of the full Hamiltonian TL. Zj is the corre- 
sponding partition function. The expectation value of an 
observable O at time t is given by 

{0)(t)=Tt{0(t)l[p j } (20) 
3 

where 0{t) is the operator associated to the observable O 
in the Heisenberg picture. 

In particular, in the following we will pay some atten- 
tion on the two-temperatures problem 

p(0) = PsPb (21) 

with a small system of size N initially prepared at a tem- 
perature T s putted in contact at one end with a big system 
of size M (called the bath) at a temperature T},. After the 
contact, the time evolution will be governed by the Hamil- 
tonian 0) with L = N + M. The expectation value takes 
the form 

(0)(t) = ^—TY{0{t)e-^ n ^e- PbUb } (22) 
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with O(t) 
is split into 



JHt 



Oe 



-VHt 



and where the total Hamiltonian 



7~C = Tt s + 7~tb + "H 1 



with 



H s = - 



l^ 1 



N 



fc=l fe=l 



^ N+M-l , N+M 

fc=iV+l fe=iV+l 



(23) 
(24) 
(25) 



and where 7i 7 = — 30^0^+1 is the initial interface in- 
teraction term. For exemple, in the situation where the 
system consists of only one spin, expressing the bath part 
in terms of its diagonal fermionic operators rf , r], one has 
for the total Hamiltonian 

h M 1 M 

u = + E £ W^ E + 1*) > ( 26 ) 

9=1 



9=1 



where 4> q (l) refers for the first bath site. The total Hamil- 
tonian thus obtained is the fermionic bath version of the 
so-called "spin-boson" model jlUlllj . The bath degrees of 
freedom act on the two-state system by the linear coupling 
-\afi = -|of EqLiKMH +Vg)- The present study 
can be interpreted in this context. 

To evaluate the expectation value (0)(t), one first has 
to express the operator 0(t) in terms of the elementary 
timc-independant r% Clifford operators. Then, one will 
have to evaluate expressions of the form 



1 



-0 b H t 



} 



(27) 



which can be calculated explicitely noticing that the 
Hamiltonians TL S and H.b can be expressed in terms of 
the lattice operators r k . 



4 Transverse magnetization 
4.1 System plus bath initial state 

The transverse magnetization is very simply expressed in 
terms of the Clifford operators, since it is a local quantity. 
One has in the Heisenberg picture 



a z n (t) = -%r^t)r^(t) 



(28) 



The expectation value of the transverse magnetization, 
with initial state p(0) — p s pb, where the first N spins of 
the chain are prepared in the canonical state p s while the 
M remaining spins are in pb, is 



2 N+M 

■ E E ( r kK(t))(n:K(t)) 

i/,v' = l k.k' — 1 



xTr{p aPb rZr£,} . (29) 



The time contractions (r%\rfi (t)) are known and given 
in terms of the (f>, tp which are obtained by the diagonal- 
ization of the T matrix. What remains to be done is to 
evaluate the trace Ti{p s pbr^r^, }. One distinguishes three 
different situations: F k and Fy are in the system part, 
r k and Fy are in the bath part, and finally one is in the 
system while the other is in the bath. In the first case, 
noticing that 



r k 2 = r sk Vfc = l...JV 



(30) 



where the r s>k s are the Clifford operators associated to 
the first part of the chain (disconnected from the rest) of 
size N, one has for the trace term 

Tr{p s p b r k rtf } = Tr s { Ps r» k r^ k , } (31) 

The same considerations, together with the fact that 

(32) 



1,2 
k+N 



where the -T^fcS are associated to the bath part and where 
Q s — °j) 1S the charge operator of the system 

part, leads to the expression 

Ti{p sPb r k+N r^ +N } = Trb{pbr^ k r( k ,} (33) 

with the indices k and k' running from 1 to M. Finally, in 
the last case with one Clifford operator in each sector, it is 
easy to show that the trace is vanishing since it is propor- 
tional to terms of the form Trh{e~ /3l,l_ ^ i7! '" 7 ^«7^ (J } = 0. 
Taking into account all the various contributions, one ob- 
tains for the expectation value of the transverse magneti- 
zation 

(**)(t)=S n (t)+B n (t) (34) 

where 



s n (t)= ]T [(n\rl(t)){ry\r^t)) 

k,k' = l 

-{r k \r^t))(ry\r^t))](s) k , k , (35) 



and 



B n (t)= Y, m + N\rm(ri, +N \rZ(t)} 



k,k' = l 



- (n +N \r^t)}(r^ +N K(t))} (B) M , (36) 

with the canonical contractions of the system 

N 

(S) M ; = -^^ 9 (fc)^ 9 (fc')tanh(/3 s -^) (37) 
9=1 

and the canonical contractions of the bath 

M 

(B) fc , fc , =-J2 h,M^ q {k') tanh (fr2f) ■ (38) 
9=1 



4 
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4.2 General expressions 

The matrices S and B contain the initial state properties. 
They can be rewritten more explicitely in terms of the 
spin operators, utilising eq. J5J) and eq. as 



(S) 



fc.fc' 



Tr{ Pa ir} k rl k ,} = -Tr{p s i(4Q%Q s kl al,}(39) 



and 



-Tr{p b ia* N+k Q b k Ql,al +k ,} 



(40) 



fc'c-^andQ^n-r 1 

the (N + M) x (N + M) matrix I as 



with Q% = nP(-^l) and Ql = Y\ K J= ^ N+] ). Defining 



S 
B 



(41) 



and the time-dependent element 

= (n\rto)){i% \n,{t))-{rl\rl{t)){rl,\r^t)) , 

it is possible to write the expectation value of the trans- 
verse magnetization as 



L—N+M 

«)(t)= E ^'WWW 
fc,fc'=l 



(43) 



It is easy to generalize that expression for initial prod- 



uct canonical states of the form 11^=1 Pj since in this case 
the initial matrix I takes the form 

/h ... 0\ 
' h 



Vo 



'■. 

I K J 



(44) 



where Ij is the initial matrix of the ^-subsystem. The 
transverse magnetization is given by (|43[1 . 
On the same lines, we have 



where all the spins are thermalized independently with 
local Hamiltonian 7i n = —\&ni the initial matrix I is 
diagonal: 

Ifc, fc < =6 k , k ,(a%){0) (48) 

where (cf)(0) is the initial transverse magnetization. In 
this case, the transverse magnetization at later time is 
given by 



In the thermodynamical limit L — > oo, the contractions 
(r k \r%(t)) depend only on the difference n — k so that the 
transverse magnetization is given by a discrete convolution 
product 



«>(*) 



k 



F t (n-k)(a* k )(0) 



(50) 



where the Green function F t (n — k) = P*'*. At the critical 
point, h = 1, the basic contractions are expressed in terms 
of Bessel functions and one can show that |12ll| 



Ft(p) = Jl p (2t) - J 2p +i(2i)J 2p _ 1 (2t) 



(51) 



In the continuum limit, one can show that the Green func- 
tion F t (p) — t/(?) with the scaling function ^2] 



/(«) 



u < 1 



(52) 



and zero otherwise. Depending on the initial state, one 
will have different profiles for the transverse magnetiza- 
tion which can be in principle explicitely obtained by per- 
forming the above convolution product. It is interesting to 
notice that the time behaviour of the transverse magneti- 
zation in this case is basically the same as in the situation 
when the starting initial state is a pure state of the form 



\W) = | . . . CTfeCTfe + l . . .) 



(53) 



where a k is the value of the z-component spin at site k 
|9I12| . One can see Ref. |1I12| for more details and physical 
examples. 



k,k'=l 



and 



'n u n+ 



i>(*)= E Ciwm 



k,k> 



fe,fe'=l 



for the bond correlations. 

5 Transverse magnetization profiles 
5.1 Completely factorized initial state 

In the completely factorized initial state 

p(0) = p? l 



(45) 



(46) 



(47) 



5.2 Two-temperature state 

We turn now to the more complex case where part of the 
chain is thermalized at inverse temperature f3 s and the rest 
at inverse temperature The transverse magnetization 
is given by the general formula (|43(l with an initial value 
matrix I which has non-diagonal terms, taking into ac- 
count that within the system part, the spins are interact- 
ing with each other and consequently have non-vanishing 
correlations, the same beeing true for the bath-part of the 
chain. 

In the following, we consider two different situations. 
The first is concerned with the interface behaviour, that is 
mainly how the jump of the transverse magnetization at 
the boundary between the system and the bath parts will 
spread out. The second situation deals with the case of a 
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finite system, of N interacting spins, initially at low tem- 
perature T s in contact at both ends with infinite chains 
at high temperature Xt>. In this case, the initially large 
transverse magnetization of the system (see the Hamilto- 
nian Q) will decay toward the small bath value as time 
is evolving. 

Consider the first situation. In order to avoid boundary 
terms, we will take the thermodynamical limit where, both 
system and bath sizes are sent to infinity. The expectation 
value of the transverse magnetization is given by 

«>(*) = x;f?(»-fc)i M -K (54) 

i,k 

where at the critical field value h = 1 , the Green functions 
are expressed in terms of Bessel functions: 

Flip) = (-If [J 2p (2t)j 2(p _ t) (2t) - J 2p+1 (2t)J 2(p _ i) _ 1 (2t)] 

The relation giving the transverse magnetization decom- 
poses into a sum of discrete convolution products 

«)(*) = 52(Ft*#)(n) (56) 

i 

with Green functions F'l defined above. In the extreme 
case where T s — and T& = oo, we arrive at a further sim- 
plification since in the initial I matrix (|4 1 11 - the B matrix 
is vanishing and the elements of the system matrix S are 
simply (S)ij = - <f>s,q(i)ips,gU)- In figure 1 we show 

the rescaled transverse magnetization profile obtained nu- 
merically by exact diagonalisation of the Hamiltonians Ji s , 
Hb and TL. 

0.75 | . 1 . 1 1 1 1 1 . 1 . 1 



0.5 - 

0.25 - 

- 

-0.25 - 

-0.5 - 




-1.5 -1 -0.5 0.5 1 1.5 

n/t 

Fig. 1. Scaling function of the transverse magnetization ob- 
tained for times t = 30, 60, 90, 120 on chains of total size 
L = 200. 



We have the very simple expression: 



where the scaling function f(u) is given by 

fi u<-l 
f{u) = < -hi -\<u<\ (58) 

1-3 «>* 

One may notice that since the excitations travel with ve- 
locity c(h = 1) = 1, one expects effects only in the causal 
region — 1 < £/t < 1 where £ is the distance from the ini- 
tial interface. Outside the causal region, both system and 
bath parts are stil equilibrated and nothing is changed 
since they are at equilibrium. The remarkable linearity of 
the profile is lost as soon as we have a departure from in- 
finite and vanishing temperature. For finite temperature 
cases, the profile is rounded. 

In the off-critical region, h > 1, the transverse mag- 
netization profile, in the case T s = and T\, = oo, shows 
a feature that was already remarked in the XX-chain [T7] 
context. Indeed, as one can see on figure 2 the magnetiza- 
tion relaxes in quantized steps. That is, if one concentrates 
near the front entering into the bath part (which has a 
vanishing magnetization, since it is at infinite tempera- 
ture), one finds a staircase like structure with constant 
area steps. If one defines 

m(£,t) = (a!)(t) - (a z bath ) , (59) 

as it is clearly shown in figure 2, the magnetization devia- 
tion from the bath magnetization {(J^ ath ) = has a scaling 
behaviour of the form 

m(l,t)=t- l '*g(!0\ (60) 

where again £ measures the distance from the initial in- 
terface, that is the interface between the system and the 
bath. Numerically, as seen from figure (J2J, it seems that 
the enveloppe of the scaling function g is a simple square- 
root, so that we have the conjecture 

m(£,t) = At-V 3 J^- £<t, (61) 

where A is a function of the field strength h. 

The scaling form l|6U|) implies that since the width of 
the steps increases as t 1 ^ 3 and their height decreases as 
t~ 1/3 , the area of the steps is indeed constant during time. 
Each step carries a definite magnetic moment (J>(h) de- 
pending on the field strength h. As h is increased toward 
infinity, the magnetic moment carried by each step reaches 
the value one (found numerically) , which is reminiscent of 
the fact that in this case, an up spin from the system ini- 
tial state \<F) — | f| ... 1") (see the Hamiltonian (JIJ) is 
injected into the bath and carried by one step. 

This simple particle-like picture does not apply at the 
critical field value h = 1, since in this case, the initial 
system state is a critical one with long range correlations. 
We have seen numerically that the quantized steps picture 
is not present in the case h < 1. 

Let us turn now to the finite system of size N, initially 
at vanishing temperature, in contact at both ends with in- 
finite chains at infinite temperature. We will refer to this 
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-1 1 — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — 1 

-14 -12 -10 -8 -6,„ -4-2 2 

(l-W 

Fig. 2. Scaling of the transverse magnetization deviation 
m(£, i) obtained for times t = 80, 120, 160, 200 at transverse 
field h = 10. 



middle of the system. The moving front between the flat 
finite transverse magnetization region and the zero mag- 
netization region is linear with a slope ~ 1/ (3t) and has 
a constant width of size N. This leads, in the asymptotic 
regime t ^> N, to the scaling behaviour 

m(n,t) ee «}(*) = M tot ij7 (^) , (62) 

where n is the site label measured from the middle of 
the system and where M to t — J2 n m ( n , t) is the total con- 
served transverse magnetization M to t — Y] n m(n, t 3> N). 
Numerically we found evidences for M tot = Nml^, v k s , 
where m^,^ = 2/tt is the T = equilibrium system bulk 
transverse magnetization. 

Finally in figure 4, we show the results obtained in the 
droplet problem in the case where the bath part is at finite 
temperature. We have more or less the same scenario as in 



situation as the droplet situation. In figure (3), we show 
the rescaled transverse magnetization profile obtained nu- 
merically at the critical field value h = 1. We see two dif- 



15 



10 - 



5 - 



















\ \\\ \\ 






n/t 



Fig. 4. Rescaled transverse magnetization with T s — 0, Tt — 1 
and at ft = 1 for different times t. Inset: Transverse magneti- 
zation profiles at different times. 



Fig. 3. Rescaled transverse magnetization with T 3 — 0, 
Tt = oo and at h = 1 for different times t. Inset: Transverse 
magnetization profiles at different times. 



ferent regimes. In the first stage, for times smaller than the 
system size N, at each interface we have the phenomenol- 
ogy depicted previously in the infinite system case, that is 
a linear behaviour of the interface given by equations l|57f) 
and l|58(l . One may also remark that some transverse mag- 
netization is lost during the initial time t ~ N . Indeed, 
since the system at t — is coupled with the bath chains, 
and since the total energy is conserved, a portion of the 
transverse energy term, e z — —h/2a z , will relax into the 
coupling energy term, e xx = —l/2a x a x : leading to a net 
loss of transverse magnetization. For times larger than the 
typical time r = TV, no more transverse magnetization is 
lost and it behaves like a conserved quantity. From the nu- 
merical profiles we see the emergence of a flat profile in the 



the infinite temperature bath case, that is a velocity one 
propagating front, which here is not linear, of constant 
size N . In the large time limit, t ^> N , numerically we 
obtain the scaling behaviour 



z,bath i ~ ' tt 

m{n,t)=m b ' + —77 



N 



(63) 



with a temperature dependent function k(u). Again, the 
total transverse magnetization is conserved in the asymp- 
totic limit t ^> N for the same reasons as discussed 
previously. If one defines the excess total magnetization 
as AMtot = SrJ m ( n >* ~> N) — ml' bath ], one obtains 
by integration of the profile l|63|l the relation AM to t = 
Y J . k(u)du. In order to recover the scaling form l|62|) . 
one has, in the limit Tf, — > oo, k(u) — > 1. In figure 5 we 
show the scaling function A; (it) obtained for different tem- 
peratures of the bath. We see clearly the flatening of the 
profile as the temperature increases. 



T. Platini and D. Karevski: Scaling and front dynamics in Ising quantum chains 



7 



0.8 




-1 -0.5 0.5 1 

U 



Fig. 5. Scaling function k(u) obtained for T b = oo, T b = 5, 
T b = 2,T b = 1, T b = 1/2, from top to bottom. 

6 Discussion 

We have investigated the front propagation and scaling 
profiles of transverse magnetization inhomogeneities in the 
Ising quantum chain. The inhomogeneities were generated 
initially by local equilibration with several thermal baths 
at different temperatures. We have concentrated our study 
on two distinct temperature configurations: the first one 
beeing the interface problem, where half of the chain is 
at a given temperature, the other half at an other tem- 
perature. The second situation considered is one, referred 
as the droplet problem, where initially a small part of 
the infinite chain is at a different temperature from the 
remaining part. In both cases, in the asymptotic time 
regime the transverse magnetization exhibits a scaling 
form t^ 1 M(n/t), where the scaling function M depends on 
the initial temperatures and on the transverse field value 
h. At the critical field value, h = 1, in the extreme case 
T\ = and T2 — 00, the transverse magnetization front 
is linear, while whitin the droplet a flat profile emerges. 
Asymptotically, in the droplet situation the scaling profile 
is given by a characteristic function 77 (u). Another inter- 
esting feature of such relaxation problem is the scaling 
behaviour of the front itself. In particular, we have shown 
that in the large field situation, h > 1, the front has a 
staircase structure with scaling i _1 / 3 g((£ — t)/t 1 / 3 ), mean- 
ing that the area of each step is conserved which implies 
a quantized relaxation of the magnetization, as already 
found in the XX-chain in Ref. |17j . 
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